Nuclear exchange generates population diversity in the wheat leaf rust pathogen Puccinia triticina.

In clonally reproducing dikaryotic rust fungi, non-sexual processes such as somatic nuclear exchange are postulated to play a role in diversity but have been difficult to detect due to the lack of genome resolution between the two haploid nuclei. We examined three nuclear-phased genome assemblies of Puccinia triticina, which causes wheat leaf rust disease. We found that the most recently emerged Australian lineage was derived by nuclear exchange between two pre-existing lineages, which originated in Europe and North America. Haplotype-specific phylogenetic analysis reveals that repeated somatic exchange events have shuffled haploid nuclei between long-term clonal lineages, leading to a global P. triticina population representing different combinations of a limited number of haploid genomes. Thus, nuclear exchange seems to be the predominant mechanism generating diversity and the emergence of new strains in this otherwise clonal pathogen. Such genomics-accelerated surveillance of pathogen evolution paves the way for more accurate global disease monitoring.

Improved quantification of Fusarium pseudograminearum (Fusarium crown rot) using qPCR measurement of infection in multi-species winter cereal experiments.

Fusarium crown rot (FCR) causes significant grain yield loss in winter cereals around the world. Breeding for resistance and/or tolerance to FCR has been slow with relatively limited success. In this study, multi-species experiments were used to demonstrate an improved method to quantify FCR infection levels at plant maturity using quantitative PCR (qPCR), as well as the genotype yield retention using residual regression deviation. Using qPCR to measure FCR infection allowed a higher degree of resolution between genotypes than traditional visual stem basal browning assessments. The results were consistent across three environments with different levels of disease expression. The improved measure of FCR infection along with genotype yield retention allows for partitioning of both tolerance and partial resistance. Together these methods offer new insights into FCR partial resistance and its relative importance to tolerance in bread wheat and barley. This new approach offers a more robust, unbiased way to select for both FCR traits within breeding programs. Key message: Genetic gain for tolerance and partial resistance against Fusarium crown rot (FCR) in winter cereals has been impeded by laborious and variable visual measures of infection severity. This paper presents results of an improved method to quantify FCR infection that are strongly correlated to yield loss and reveal previously unrecognised partial resistance in barley and wheat varieties.

A thousand-genome panel retraces the global spread and adaptation of a major fungal crop pathogen.

Human activity impacts the evolutionary trajectories of many species worldwide. Global trade of agricultural goods contributes to the dispersal of pathogens reshaping their genetic makeup and providing opportunities for virulence gains. Understanding how pathogens surmount control strategies and cope with new climates is crucial to predicting the future impact of crop pathogens. Here, we address this by assembling a global thousand-genome panel of Zymoseptoria tritici, a major fungal pathogen of wheat reported in all production areas worldwide. We identify the global invasion routes and ongoing genetic exchange of the pathogen among wheat-growing regions. We find that the global expansion was accompanied by increased activity of transposable elements and weakened genomic defenses. Finally, we find significant standing variation for adaptation to new climates encountered during the global spread. Our work shows how large population genomic panels enable deep insights into the evolutionary trajectory of a major crop pathogen.

Multi-stage resistance to Zymoseptoria tritici revealed by GWAS in an Australian bread wheat diversity panel.

Septoria tritici blotch (STB) has been ranked the third most important wheat disease in the world, threatening a large area of wheat production. Although major genes play an important role in the protection against Zymoseptoria tritici infection, the lifespan of their resistance unfortunately is very short in modern wheat production systems. Combinations of quantitative resistance with minor effects, therefore, are believed to have prolonged and more durable resistance to Z. tritici. In this study, new quantitative trait loci (QTLs) were identified that are responsible for seedling-stage resistance and adult-plant stage resistance (APR). More importantly was the characterisation of a previously unidentified QTL that can provide resistance during different stages of plant growth or multi-stage resistance (MSR). At the seedling stage, we discovered a new isolate-specific QTL, QSt.wai.1A.1. At the adult-plant stage, the new QTL QStb.wai.6A.2 provided stable and consistent APR in multiple sites and years, while the QTL QStb.wai.7A.2 was highlighted to have MSR. The stacking of multiple favourable MSR alleles was found to improve resistance to Z. tritici by up to 40%.

Remarkable recent changes in the genetic diversity of the avirulence gene AvrStb6 in global populations of the wheat pathogen Zymoseptoria tritici.

Septoria tritici blotch (STB), caused by the fungus Zymoseptoria tritici, is one of the most economically important diseases of wheat. Recently, both factors of a gene-for-gene interaction between Z. tritici and wheat, the wheat receptor-like kinase Stb6 and the Z. tritici secreted effector protein AvrStb6, have been identified. Previous analyses revealed a high diversity of AvrStb6 haplotypes present in earlier Z. tritici isolate collections, with up to c.18% of analysed isolates possessing the avirulence isoform of AvrStb6 identical to that originally identified in the reference isolate IPO323. With Stb6 present in many commercial wheat cultivars globally, we aimed to assess potential changes in AvrStb6 genetic diversity and the incidence of haplotypes allowing evasion of Stb6-mediated resistance in more recent Z. tritici populations. Here we show, using targeted resequencing of AvrStb6, that this gene is universally present in field isolates sampled from major wheat-growing regions of the world in 2013-2017. However, in contrast to the data from previous AvrStb6 population studies, we report a complete absence of the originally described avirulence isoform of AvrStb6 amongst modern Z. tritici isolates. Moreover, a remarkably small number of haplotypes, each encoding AvrStb6 protein isoforms conditioning virulence on Stb6-containing wheat, were found to predominate among modern Z. tritici isolates. A single virulence isoform of AvrStb6 was found to be particularly abundant throughout the global population. These findings indicate that, despite the ability of Z. tritici to sexually reproduce on resistant hosts, AvrStb6 avirulence haplotypes tend to be eliminated in subsequent populations.

The identification of a transposon affecting the asexual reproduction of the wheat pathogen Zymoseptoria tritici.

Zymoseptoria tritici, the causal agent of Septoria tritici blotch, is a fungal wheat pathogen that causes significant global yield losses. Within Z. tritici populations, quantitative differences in virulence among different isolates are commonly observed; however, the genetic components that underpin these differences remain elusive. In this study, intraspecific comparative transcriptomic analysis was used to identify candidate genes that contribute to differences in virulence on the wheat cultivar WW2449. This led to the identification of a multicopy gene that was not expressed in the high-virulence isolate when compared to the medium- and low-virulence isolates. Further investigation suggested this gene resides in a 7.9-kb transposon. Subsequent long-read sequencing of the isolates used in the transcriptomic analysis confirmed that this gene did reside in an active Class II transposon, which is composed of four genes named REP9-1 to -4. Silencing and overexpression of REP9-1 in two distinct genetic backgrounds demonstrated that its expression alone reduces the number of pycnidia produced by Z. tritici during infection. The REP9-1 gene identified within a Class II transposon is the first discovery of a gene in a transposable element that influences the virulence of Z. tritici. This discovery adds further complexity to genetic loci that contribute to quantitative virulence in this important pathogen.

The complex genomic basis of rapid convergent adaptation to pesticides across continents in a fungal plant pathogen.

Convergent evolution leads to identical phenotypic traits in different species or populations. Convergence can be driven by standing variation allowing selection to favour identical alleles in parallel or the same mutations can arise independently. However, the molecular basis of such convergent adaptation remains often poorly resolved. Pesticide resistance in agricultural ecosystems is a hallmark of convergence in phenotypic traits. Here, we analyse the major fungal pathogen Zymoseptoria tritici causing serious losses on wheat and with fungicide resistance emergence across several continents. We sampled three population pairs each from a different continent spanning periods early and late in the application of fungicides. To identify causal loci for resistance, we combined knowledge from molecular genetics work and performed genome-wide association studies (GWAS) on a global set of isolates. We discovered yet unknown factors in azole resistance including a gene encoding membrane associated functions. We found strong support for the "hotspot" model of resistance evolution with convergent changes in a small set of loci but additional loci showed more population-specific allele frequency changes. Genome-wide scans of selection showed that half of all known resistance loci were overlapping a selective sweep region. Hence, the application of fungicides was one of the major selective agents acting on the pathogen over the past decades. Furthermore, loci identified through GWAS showed the highest overlap with selective sweep regions underlining the importance to map phenotypic trait variation in evolving populations. Our population genomic analyses highlighted that both de novo mutations and gene flow contributed to convergent pesticide adaptation.

Recent insights into barley and Rhynchosporium commune interactions.

Rhynchosporium commune is the causal pathogen of scald in barley (Hordeum vulgare), a foliar disease that can reduce yield by up to 40% in susceptible cultivars. R. commune is found worldwide in all temperate growing regions and is regarded as one of the most economically important barley pathogens. It is a polycyclic pathogen with the ability to rapidly evolve new virulent strains in response to resistance genes deployed in commercial cultivars. Hence, introgression and pyramiding of different loci for resistance (qualitative or quantitative) through marker-assisted selection is an effective way to improve scald resistance in barley. This review summarizes all 148 resistance quantitative trait loci reported at the date of submission of this review and projects them onto the barley physical map, where it is clear many loci co-locate on chromosomes 3H and 7H. We have summarized the major named resistance loci and reiterated the renaming of Rrs15 (CI8288) to Rrs17. This review provides a comprehensive resource for future discovery and breeding efforts of qualitative and quantitative scald resistance loci.

Transposon-Mediated Horizontal Transfer of the Host-Specific Virulence Protein ToxA between Three Fungal Wheat Pathogens.

Most known examples of horizontal gene transfer (HGT) between eukaryotes are ancient. These events are identified primarily using phylogenetic methods on coding regions alone. Only rarely are there examples of HGT where noncoding DNA is also reported. The gene encoding the wheat virulence protein ToxA and the surrounding 14 kb is one of these rare examples. ToxA has been horizontally transferred between three fungal wheat pathogens (Parastagonospora nodorum, Pyrenophora tritici-repentis, and Bipolaris sorokiniana) as part of a conserved ∼14 kb element which contains coding and noncoding regions. Here we used long-read sequencing to define the extent of HGT between these three fungal species. Construction of near-chromosomal-level assemblies enabled identification of terminal inverted repeats on either end of the 14 kb region, typical of a type II DNA transposon. This is the first description of ToxA with complete transposon features, which we call ToxhAT. In all three species, ToxhAT resides in a large (140-to-250 kb) transposon-rich genomic island which is absent in isolates that do not carry the gene (annotated here as toxa- ). We demonstrate that the horizontal transfer of ToxhAT between P. tritici-repentis and P. nodorum occurred as part of a large (∼80 kb) HGT which is now undergoing extensive decay. In B. sorokiniana, in contrast, ToxhAT and its resident genomic island are mobile within the genome. Together, these data provide insight into the noncoding regions that facilitate HGT between eukaryotes and into the genomic processes which mask the extent of HGT between these species.IMPORTANCE This work dissects the tripartite horizontal transfer of ToxA, a gene that has a direct negative impact on global wheat yields. Defining the extent of horizontally transferred DNA is important because it can provide clues to the mechanisms that facilitate HGT. Our analysis of ToxA and its surrounding 14 kb suggests that this gene was horizontally transferred in two independent events, with one event likely facilitated by a type II DNA transposon. These horizontal transfer events are now in various processes of decay in each species due to the repeated insertion of new transposons and subsequent rounds of targeted mutation by a fungal genome defense mechanism known as repeat induced point mutation. This work highlights the role that HGT plays in the evolution of host adaptation in eukaryotic pathogens. It also increases the growing body of evidence indicating that transposons facilitate adaptive HGT events between fungi present in similar environments and hosts.

Rapid Parallel Evolution of Azole Fungicide Resistance in Australian Populations of the Wheat Pathogen Zymoseptoria tritici.

Zymoseptoria tritici is a globally distributed fungal pathogen which causes Septoria tritici blotch on wheat. Management of the disease is attempted through the deployment of resistant wheat cultivars and the application of fungicides. However, fungicide resistance is commonly observed in Z. tritici populations, and continuous monitoring is required to detect breakdowns in fungicide efficacy. We recently reported azole-resistant isolates in Australia; however, it remained unknown whether resistance was brought into the continent through gene flow or whether resistance emerged independently. To address this question, we screened 43 isolates across five Australian locations for azole sensitivity and performed whole-genome sequencing on 58 isolates from seven locations to determine the genetic basis of resistance. Population genomic analyses showed extremely strong differentiation between the Australian population recovered after azoles began to be used and both Australian populations recovered before azoles began to be used and populations on different continents. The apparent absence of recent gene flow between Australia and other continents suggests that azole fungicide resistance has evolved de novo and subsequently spread within Tasmania. Despite the isolates being distinct at the whole-genome level, we observed combinations of nonsynonymous substitutions at the CYP51 locus identical to those observed elsewhere in the world. We observed nine previously reported nonsynonymous mutations as well as isolates that carried a combination of the previously reported L50S, S188N, A379G, I381V, Y459DEL, G460DEL, and N513K substitutions. Assays for the 50% effective concentration against a subset of isolates exposed to the tebuconazole and epoxiconazole fungicides showed high levels of azole resistance. The rapid, parallel evolution of a complex CYP51 haplotype that matches a dominant European haplotype demonstrates the enormous potential for de novo resistance emergence in pathogenic fungi.IMPORTANCE Fungicides are essential to control diseases in agriculture because many crops are highly susceptible to pathogens. However, many pathogens rapidly evolve resistance to fungicides. A large body of studies have described specific mutations conferring resistance and have often made inferences about the origins of resistance based on sequencing data from the target gene alone. Here, we show the de novo acquisition of resistance to the ubiquitously used azole fungicides in genetically isolated populations of the wheat pathogen Zymoseptoria tritici in Tasmania, Australia. We confirm evidence for parallel evolution through genome-scale analyses of representative worldwide populations. The emergence of complex resistance haplotypes following a well-documented recent introduction of azoles into Australian farming practices demonstrates how rapidly chemical resistance evolves in agricultural ecosystems.

Bivariate analysis of barley scald resistance with relative maturity reveals a new major QTL on chromosome 3H.

The disease scald of barley is caused by the pathogen Rhynchosporium commune and can cause up to 30-40% yield loss in susceptible cultivars. In this study, the Australian barley cultivar 'Yerong' was demonstrated to have resistance that differed from Turk (Rrs1 (Rh3 type)) based on seedling tests with 11 R. commune isolates. A doubled haploid population with 177 lines derived from a cross between 'Yerong' and the susceptible Australian cultivar 'Franklin' was used to identify quantitative trait loci (QTL) for scald resistance. A QTL on chromosome 3H was identified with large effect, consistent with a major gene conferring scald resistance at the seedling stage. Under field conditions, a bivariate analysis was used to model scald percentage of infected leaf area and relative maturity, the residuals from the regression were used as our phenotype for QTL analysis. This analysis identified one major QTL on chromosome 3H, which mapped to the same position as the QTL at seedling stage. The identified QTL on 3H is proposed to be different from the Rrs1 on the basis of seedling resistance against different R. commune isolates and physical map position. This study increases the current understanding of scald resistance and identifies genetic material possessing QTLs useful for the marker-assisted selection of scald resistance in barley breeding programs.

Genetic mapping of Stb19, a new resistance gene to Zymoseptoria tritici in wheat.

KEY MESSAGE: A new and dominant R gene Stb19 is identified from a soft wheat cultivar 'Lorikeet' and was mapped on the distal region of chromosome 1DS. Two tightly linked KASP markers were also discovered and validated for molecular-assisted breeding programs. A new R gene, designated as Stb19, provides resistance to Zymoseptoria tritici in wheat. This new dominant gene resides on the short arm of chromosome 1D, exhibiting complete resistance to three Z. tritici isolates, WAI332, WAI251, and WAI161, at the seedling stage. A genetic linkage map, based on an F2:3 population of 'Lorikeet' and 'Summit,' found the Stb19 gene at a 9.3 cM region on 1DS, closely linked with two Kompetitive Allele-Specific PCR markers, snp_4909967 and snp_1218021. Further, the two markers were tested and validated in another F2:3 population and 266 different wheat accessions, which gave over 95% accuracy of resistance/susceptibility prediction. Combined with the physical location of the identified SNPs and the previous evidence of gene order on chromosome 1DS (centromere-Sr45-Sr33-Lr21-telomere), Stb19 is proposed to be located between Sr33 and Lr21. Thus, the newly discovered Stb19 along with the KASP markers represents an increase in genetic resources available for wheat breeding resistance to Z. tritici.

Accounting for Genotype-by-Environment Interactions and Residual Genetic Variation in Genomic Selection for Water-Soluble Carbohydrate Concentration in Wheat.

Abiotic stress tolerance traits are often complex and recalcitrant targets for conventional breeding improvement in many crop species. This study evaluated the potential of genomic selection to predict water-soluble carbohydrate concentration (WSCC), an important drought tolerance trait, in wheat under field conditions. A panel of 358 varieties and breeding lines constrained for maturity was evaluated under rainfed and irrigated treatments across two locations and two years. Whole-genome marker profiles and factor analytic mixed models were used to generate genomic estimated breeding values (GEBVs) for specific environments and environment groups. Additive genetic variance was smaller than residual genetic variance for WSCC, such that genotypic values were dominated by residual genetic effects rather than additive breeding values. As a result, GEBVs were not accurate predictors of genotypic values of the extant lines, but GEBVs should be reliable selection criteria to choose parents for intermating to produce new populations. The accuracy of GEBVs for untested lines was sufficient to increase predicted genetic gain from genomic selection per unit time compared to phenotypic selection if the breeding cycle is reduced by half by the use of GEBVs in off-season generations. Further, genomic prediction accuracy depended on having phenotypic data from environments with strong correlations with target production environments to build prediction models. By combining high-density marker genotypes, stress-managed field evaluations, and mixed models that model simultaneously covariances among genotypes and covariances of complex trait performance between pairs of environments, we were able to train models with good accuracy to facilitate genetic gain from genomic selection.

The discovery of the virulence gene ToxA in the wheat and barley pathogen Bipolaris sorokiniana.

Bipolaris sorokiniana is the causal agent of multiple diseases on wheat and barley and is the primary constraint to cereal production throughout South Asia. Despite its significance, the molecular basis of disease is poorly understood. To address this, the genomes of three Australian isolates of B. sorokiniana were sequenced and screened for known pathogenicity genes. Sequence analysis revealed that the isolate BRIP10943 harboured the ToxA gene, which has been associated previously with disease in the wheat pathogens Parastagonospora nodorum and Pyrenophora tritici-repentis. Analysis of the regions flanking ToxA within B. sorokiniana revealed that it was embedded within a 12-kb genomic element nearly identical to the corresponding regions in P. nodorum and P. tritici-repentis. A screen of 35 Australian B. sorokiniana isolates confirmed that ToxA was present in 12 isolates. Sequencing of the ToxA genes within these isolates revealed two haplotypes, which differed by a single non-synonymous nucleotide substitution. Pathogenicity assays showed that a B. sorokiniana isolate harbouring ToxA was more virulent on wheat lines that contained the sensitivity gene when compared with a non-ToxA isolate. This work demonstrates that proteins that confer host-specific virulence can be horizontally acquired across multiple species. This acquisition can dramatically increase the virulence of pathogenic strains on susceptible cultivars, which, in an agricultural setting, can have devastating economic and social impacts.

Selection for water-soluble carbohydrate accumulation and investigation of genetic × environment interactions in an elite wheat breeding population.

KEY MESSAGE: Water-soluble carbohydrate accumulation can be selected in wheat breeding programs with consideration of genetic × environmental interactions and relationships with other important characteristics such as relative maturity and nitrogen concentration, although the correlation between WSC traits and grain yield is low and inconsistent. The potential to increase the genetic capacity for water-soluble carbohydrate (WSC) accumulation is an opportunity to improve the drought tolerance capability of rainfed wheat varieties, particularly in environments where terminal drought is a significant constraint to wheat production. A population of elite breeding germplasm was characterized to investigate the potential for selection of improved WSC concentration and total amount in water deficit and well-watered environments. Accumulation of WSC involves complex interactions with other traits and the environment. For both WSC concentration (WSCC) and total WSC per area (WSCA), strong genotype × environment interactions were reflected in the clear grouping of experiments into well-watered and water deficit environment clusters. Genetic correlations between experiments were high within clusters. Heritability for WSCC was larger than for WSCA, and significant associations were observed in both well-watered and water deficit experiment clusters between the WSC traits and nitrogen concentration, tillering, grains per m2, and grain size. However, correlations between grain yield and WSCC or WSCA were weak and variable, suggesting that selection for these traits is not a better strategy for improving yield under drought than direct selection for yield.

Genome-Wide Associations for Water-Soluble Carbohydrate Concentration and Relative Maturity in Wheat Using SNP and DArT Marker Arrays.

Improving water-use efficiency by incorporating drought avoidance traits into new wheat varieties is an important objective for wheat breeding in water-limited environments. This study uses genome wide association studies (GWAS) to identify candidate loci for water-soluble carbohydrate accumulation-an important drought-avoidance characteristic in wheat. Phenotypes from a multi-environment trial with experiments differing in water availability and separate single nucleotide polymorphism (SNP) and diversity arrays technology (DArT) marker sets were used to perform the analyses. Significant associations for water-soluble carbohydrate accumulation were identified on chromosomes 1A, 1B, 1D, 2D, and 4A. Notably, these loci did not collocate with the major loci identified for relative maturity. Loci on chromosome 1D collocated with markers previously associated with the high molecular weight glutenin Glu-D1 locus. Genetic × environmental interactions impacted the results strongly, with significant associations for carbohydrate accumulation identified only in the water-deficit experiments. The markers associated with carbohydrate accumulation may be useful for marker-assisted selection of drought tolerance in wheat.

Utilizing Gene Tree Variation to Identify Candidate Effector Genes in Zymoseptoria tritici.

Zymoseptoria tritici is a host-specific, necrotrophic pathogen of wheat. Infection by Z. tritici is characterized by its extended latent period, which typically lasts 2 wks, and is followed by extensive host cell death, and rapid proliferation of fungal biomass. This work characterizes the level of genomic variation in 13 isolates, for which we have measured virulence on 11 wheat cultivars with differential resistance genes. Between the reference isolate, IPO323, and the 13 Australian isolates we identified over 800,000 single nucleotide polymorphisms, of which ∼10% had an effect on the coding regions of the genome. Furthermore, we identified over 1700 probable presence/absence polymorphisms in genes across the Australian isolates using de novo assembly. Finally, we developed a gene tree sorting method that quickly identifies groups of isolates within a single gene alignment whose sequence haplotypes correspond with virulence scores on a single wheat cultivar. Using this method, we have identified < 100 candidate effector genes whose gene sequence correlates with virulence toward a wheat cultivar carrying a major resistance gene.

Occurrence of Winter Cereal Viruses in New South Wales, Australia, 2006 to 2014.

Winter cereal viruses can cause significant crop losses; however, detailed knowledge of their occurrence in New South Wales, Australia is very limited. This paper reports on the occurrence of Wheat streak mosaic virus (WSMV), Wheat mosaic virus (WMoV), Barley yellow dwarf virus (BYDV), Cereal yellow dwarf virus (CYDV), and their serotypes between 2006 and 2014. Detection of WMoV is confirmed in eastern Australia for the first time. The BYDV and CYDV 2014 epidemic is examined in detail using 139 samples of wheat, barley, and oat surveyed from southern New South Wales. The presence of virus was determined using enzyme-linked immunosorbent assays. The results reveal a high frequency of the serotype Barley yellow dwarf virus - MAV as a single infection present in 27% of samples relative to Barley yellow dwarf virus - PAV in 19% and CYDV in 14%. Clear differences emerged in the infection of different winter cereal species by serotypes of BYDV and CYDV. These results are contrasted to other Australian and international studies.

Next-generation re-sequencing as a tool for rapid bioinformatic screening of presence and absence of genes and accessory chromosomes across isolates of Zymoseptoria tritici.

The wheat pathogen Zymoseptoria tritici possesses a large number of accessory chromosomes that may be present or absent in its genome. The genome of the reference isolate IPO323 has been assembled to a very high standard and contains 21 full length chromosome sequences, 8 of which represent accessory chromosomes. The IPO323 reference, when combined with low-cost next-generation sequencing and bioinformatics, can be used as a powerful tool to assess the presence or absence of accessory chromosomes. We present an outline of a range of bioinformatics techniques that can be applied to the analysis of presence-absence variation among accessory chromosomes across 13 novel isolates of Z. tritici.